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^ . Abstract 

^ I We calculate the gluon structure function of a color dipole in near-light-cone SU(2) lattice QCD 
> ■ 

ijr^ ' as a function of • The quark and antiquark are external non-dynamical degrees of freedom which 

I ! 

■ act as sources of the gluon string configuration defining the dipole. We compute the color dipole 
p . 

, matrix element of transversal chromo-electric and chromo-magnetic field operators separated along 

O 

Q\ I a direction close to the light cone, the Fourier transform of which is the gluon structure function. 



As vacuum state in the pure glue sector, we use a variational ground state of the near-light-cone 



^ , Hamiltonian. We derive a recursion relation for the gluon structure function on the lattice similar 

H ; 

_Cu ^ to the perturbative DGLAP equation. It depends on the number of transversal links assembling 

the Schwinger string of the dipole. Fixing the mean momentum fraction of the gluons to the 
"experimental value" in a proton, we compare our gluon structure function for a dipole state with 
four links with the NLO MRST 2002 and the CTEQAB-0 parameterizations at = 1.5 GeV^. 
Within the systematic uncertainty we find rather good agreement. We also discuss the low xb 
behavior of the gluon structure function in our model calculation. 
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I. INTRODUCTION 



First support for QCD as the theory of strong interactions has come from deep inelastic 
scattering. The structure of the proton unfolds itself in terms of partons which interact 
only weakly due to asymptotic freedom. The deduced structure functions representing the 
proton constituents are Fourier transforms of quark/gluon operators separated by lie 
distances. Theoretical calculations of scaling violations by the DGLAP equation [l|, 
have strongly contributed to the understanding of deep inelastic scattering. A current QCD 
analysis of experiments is given e.g. in Refs. {4,5,^. For small xb, there appear - in addition 



it like 

,33 



3,y,3 



to terms from DGLAP-evolution - other BFKL-contributions like ^^(Q^) log(l/xB) 
which need special care. Also, it is needless to stress that in a perturbative framework 
structure functions themselves cannot be calculated from first principles. Euclidean lattice 
simulations use the operator product expansion to get information about quark structure 
functions. In this way, the lowest moments of the meson and nucleon structure functions 
have been evaluated jiol . [ill . 12|. Nowadays, this method has been generalized to non- 
forward matrix elements (generalized parton distribution functions) 13j. Recently, loop- 
loop correlation functions of tilted Wegner- Wilson loops have been cornputed on a Euclidean 
lattice [14 1 which can be related to the gluon distribution function [l5, 16, 17 1 of a color 
dipole or a hadron. 

Parallel to these investigations, the light cone lattice community has pursued 



18 



the idea of a different formulation of QCD on or near the light cone. The hope is that 
a theoretical framework based on constituents moving along the light cone will be simple, 
rather closely following the course of the experimental discovery of quarks. Of course, the 
light cone approach has to attempt to incorporate the non-perturbative QCD vacuum, which 
is hard to achieve in the framework of a Fock representation of free fields acting on a trivial 
vacuum. Also rewriting a spatially quantised theory into a theory quantised on a light like 
surface may present problems. 



Therefore, in a recent paper 211] we have developed a near-light-cone (nlc) approach in 
which we can combine the advantages of the lattice formulation with the advantages of light 
cone simplifications. In this reference, we have constructed a ground state wave functional 
of the near-light-cone Hamiltonian which, in the light cone limit, becomes simpler than 
the equal time ground state in the similar strong coupling approximation. Here, in a first 
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application we use this variationally optimized ground state wave functional to determine 
the gluon distribution function of a color dipole. A color dipole is a system consisting of a 
static quark and antiquark pair connected by a Schwinger string. In our simplified dipole 
picture, we use the average momentum fraction carried by the gluons as extracted from 
phenomenological analyses as an input. We then predict the shape of the gluon structure 
function as a function of the transverse size of the dipole. One should remark that most 
dipole calculations for the gluon distribution are done in a reference system where the hadron 
under consideration is at rest. In our calculation the hadron is attached to the fast moving 
frame. This justifies the application of near-light-cone dynamics for its constituents. 

The outline of the paper is as follows: In Sec. [ITl we review the original definition of the 
gluon structure function on the light cone and some of its properties. Next, we define the 
gluon structure function near the light cone such that its original definition is recovered in 
the light cone limit. In Sec. IIIH we recapitulate properties of the near-light-cone (nlc) lattice 
formulation. This leads us to the variationally optimized ground state wave functional. In 
Sec. IIVI the model of the hadron as a color dipole state is outlined. We define the lattice 
counterpart of the nlc correlation function in Sec. |Vl In Sec. I VII we discuss the lattice 
computation of the gluon structure function for a one-link dipole, which yields the building 
block for the computation of hadronic gluon distribution functions. Sec. IVIII contains our 
results and their interpretations. In Sec. IVIIII we formulate our conclusions and discuss 
possible improvements. 

II. GLUON DISTRIBUTION FUNCTION ON AND CLOSE TO THE LIGHT 
CONE 

First, we review the original definition and some properties of the gluon distribution 
function. In deep inelastic scattering the hadronic target is probed on the light cone, i.e. 
at equal light cone time = 0. Here, = (x° + x^)/\/2 and x~ = {x'^ — x^)/^/2 
are the ordinary light cone temporal and the light cone longitudinal coordinate. In light 
cone quantisation one quantises on a light like hypersurface defined by = 0. On this 
hypersurface, the entire scattering process is static. Therefore, there is no need to evolve 
the hadronic wave function in light cone time during scattering. 

The internal structure of the hadronic target is encoded in parton distribution functions. 
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For example, the gluon distribution function g{xB) represents the probabihty that a gluon 
carries the longitudinal momentum fraction xb of the fast moving hadronic target i22| . In 
light cone coordinates, it is given by the Fourier transform of the matrix element of the two- 
point operator G{z^ , zj_ ; 0, z±) of longitudinally separated gluon field strength operators in 
a hadron state \h{p_,6j_)): 

1 1 /"^ . .o^ „ 1 



g[XB) 



2 7r 



dz- d'z^e'''-'^P-' 



P- 



/i(p_,0_l) G{z ,z±;0,z^) /i(p_,0_L 



The notation \h{p^,0±)) emphasizes that the hadron is localized with its center of mass 
in transversal configuration space at b± = Oj_ and carries longitudinal momentum p_. The 
momentum is normalized relative to the total momentum p_ of the entire hadron. The 
index "c" indicates that the connected matrix element is taken, i.e. the product of the 
vacuum matrix element with the normalization of the hadronic state is subtracted: 

(/i(p_,Ox)| G{z-,z^;0,z^) |%_,0x))^ 

= (%-,0^)| 0,il) |%-,0x)) 



n 



G{z-,z^;0,z^) n) (hip_,6^) h{p^,6^ 



h{p_,0^) h{p_,0^_)) = 2p_L 



(2) 



Here, \Q) denotes the vacuum state. L_ is the spatial extension along the longitudinal 
direction of the normalization box. The point split operator G{z^ , z± ; 0, z±) corresponding 
to the wanted correlation function is given by [22| 

2 

G{z-,z^;0,z^) = 0,zl)F!fc(0,il) (3) 



k=l 



with 



Vexp-\ig / dx A1{x , z±) X 



c 

adj 



ab 



(4) 



The Schwinger string in the adjoint representation ^^(z", z± ; 0, zj_) connects the gluon field 
strength operator at the point zj_ in the z^ = plane, F^f^{0,z±), with the longitudinally 
separated field strength operator F'^i^{z^ , z±) along a light like path. In the usual light cone 
quantisation approach, one uses the so called light cone gauge A_ = = for quanti- 
sation. This sets the Schwinger string along the light cone equal to one. The importance 
of the Schwinger strings along the light cone is visualized e.g. in the loop-loop correlation 



4 



model where hadron-hadron scattering cross-sections are calculated from Wegner- Wilson 



loop correlation functions 
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24| . The eikonal phases arising from the strings along the 



x~ direction also describe the so-called "final state" interaction effects which distinguish 
structure functions from parton probabilities 

The gluon distribution function defined in that way obeys a momentum sum rule, i.e. 
the average momentum fraction of the hadron carried by the gluons is related to the first 
moment of the structure function. The integral over xb can be reformulated as an integral 
over gluon momenta p^, which yields the matrix element of the gluonic two-point operator 
G taken at z~ = 0: 
"1 



{xb) 





1 

2^ 



dXB XBg{xB) 



P- 



dz d^z^—^5{z ) (h{p^,0) G{z ,z±;0,z±) h{p^,0±) 
2p_ \ 



(5) 



which coincides with the longitudinal light cone momentum density operator V 



Ic 



G{z ,z_i_]0, z±) 



(6) 



k=l 



Hence, the average gluon fractional momentum is given by 

1 



{xb) 



2 pi 



d'z^ (hip^,0^) P!f(0,zx) %_,0j 



(7) 



The subscript c could be dropped since the disconnected part of the matrix element vanishes. 
This is the case because the expectation value of the longitudinal momentum density of the 
vacuum state vanishes, see Sec. |Vl Eq. (I7l) yields the average fractional gluon momentum 
with the normalization of the hadronic state given in Eq. ([2]). 

To compute the gluon distribution function non-perturbatively on the lattice, we shall 
use the light cone limit of "near to the lig ht cone" (nlc) quantisation instead. Here, nlc 



refers to near-light-cone coordinates 



26 



271] which have been introduced to implement light 



front quantisation as a limit of equal time quantisation. The nlc transverse and longitudinal 
coordinates x± and x~ are defined in a similar way as usual light cone coordinates. The def- 
inition of the temporal nlc coordinate x~^ however contains an additional external parameter 
7] which parameterizes a rotation in the x° — x^ plane not included in the Lorentz group and 
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which allows for a smooth interpolation between equal time quantisation {i] = 1 , = 
and light cone quantisation [r] = , = 1/2 {x^ + x^)). 

1 



X 



X 



(1 + r^^) x° + (1 - r^^) 



x' 



x'-x^ 



(8) 



Note that the 77 ^ limit can be interpreted as the infinite momentum frame limit in which 
the partons of the color dipole move with infinite momentum. Quantisation in terms of near- 
light-cone coordinates in contrast to ordinary light cone quantisation has the advantage that 
no quantum constraint equations have to be solved. This makes a lattice treatment feasible, 
at least in the Hamiltonian formulation. 

In a Hamiltonian nlc theory obeying the standard = gauge, the operator of the 
longitudinal momentum density in the pure gauge sector can be obtained from the energy 
momentum tensor by expressing the temporal components of the field strength tensor in 
terms of the chromo-electric field operators similar to the usual Legendre transformation 
from the Yang-Mills Lagrange density to the Hamiltonian density. In the = gauge, 
the chromo-electric field operators are given by the functional derivatives of the Lagrange 
density with respect to the space-time components of the field strength tensor. Hence, the 



operator of the longitudinal momentum density is given by 2l|] 



H 2 

^-(^^) = 2 5Z (n^(^1 F-ki^ + F-ki^) m^)) ■ (9) 

k=l 

Here, H"(z) is the chromo-electric field operator which is canonically conjugate to the gauge 
field A'^jiy), i.e. 

[lin^,A%y)]=-i6^'\z-y)6,,6^'' , t,j = l,2,-. (10) 

The longitudinal momentum density in Eq. iQ is symmetrized in order to render it 
hermitean. This combination of transversal chromo-electric field operators and chromo- 
magnetic field operators is quite natural because it resembles the Poynting vector in ordinary 
electrodynamics representing the momentum density of the electromagnetic field in longi- 
tudinal direction. In the light cone limit, the transversal chro mo- magnetic field strength 
operators become equal to the corresponding chromo-electric field strength operators due to 
the constraint equation which emerges in the light cone limit. 
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In order to have the same momentum sum rule in near-hght-cone coordinates as one 
has in hght cone coordinates, we define the operator corresponding to the near-hght-cone 
correlation function as a point-split generalization of the longitudinal momentum density 
given in Eq. 

G{z-,zr,z'-,z^) = - — J2{^ti^-,^±)Sttiz-,z^;z'-z^)F\{z'~,z^) 

k=l 

+ Utiz'-, il) St,iz'-, ; z' z^) F\iz-, z^) + h.c) .(11) 



We have symmetrized also this operator with respect to an interchange z~ ^ z'^ and with 
respect to the ordering of the transversal chromo-magnetic and -electric field operators. In 
the light cone limit, this operator reproduces the definition Eq. ([T]). Note that in our nlc 
Hamiltonian approach the gauge fields A- are fully dynamical gauge fields, only the gauge 
choice = has been implemented. 

For later convenience, we use translation invariance of the expectation value in order to 
introduce an additional integration over the longitudinal coordinate z'~ . Using the above 
operator Eq. flTTl) and the normalization of the hadronic state Eq. ([2]), the gluon distribution 
function in nlc coordinates is given by 

^(xB) = lim / dz~ d:^~ d^z^e-'^^"^-^'" -^"^g{z- ,z^-z!- ,zA (12) 

r?^0 2tXXb J-oo 

with the abbreviation 





G{z ,zx_\z' ,£l) 


/i(p_,Ox)) 

/ c 




/i(p_,Ox)) 



~g{z-,z^- zl) = 2 ^ -. ^ . (13) 



Note that the longitudinal nlc momentum p_ of the target state is restricted to positive 
values in the light cone limit. If p_ is expressed in terms of the ordinary Minkowski space 
momentum ps, one obtains for an on-shell particle like the target hadron the following 
expression for the longitudinal momentum in the nlc frame: 

in n+^lb7P ^ limp_G[0,oo]. (14) 



Here, = + vn? is the transversal mass. 
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III. NEAR-LIGHT-CONE LATTICE HAMILTONIAN 



In our previous work [2l| we have regularized SU{2) gauge theory for our purposes by 
introducing a spatial nlc lattice. The size is L_L']_ = N_a_N\a\, where N_,N±_ are the 
number of lattice sites along the light like direction x~ and the transversal directions Xk- The 
lattice spacings in these directions are a_ and a±. In the following, all spatial and momentum 
variables are assumed to be made dimensionless lattice quantities by multiplication with the 
appropriate powers of the lattice spacings in longitudinal and transversal directions. The 
gauge field degrees of freedom on the lattice are given by the gluon link operators Uj{x) 



Uj (x) = V exp {^9 j + y^j] 



(15) 



The ordering symbol V orders the products of gluon fields Aj{x + y'ej) from left {y = 0) to 
right y = 1). In Hamiltonian theory we have two transversal gauge fields Aj, {j = 1, 2) and 
one longitudinal gauge field A_, (j = — ). The A+ component of the gauge field is set equal 
to zero in the Hamiltonian approach. As a result one has the Gauss-law constraint which 
restricts the entire Hilbert space to the physical sector of gauge invariant states. The gluon 
dynamics is determined by the effective nlc lattice Hamiltonian which has been derived in 
Ref. |2l|. It represents the gluon energy density on the lattice. The QCD coupling constant 
enters as A = 4/(7^ in the SU{2) Hamiltonian 



K 



efljat 




1 - Re Ui2(x 



f Im f/_fe(x) 



(16) 



The 2x2 Pauli matrices (j"/2 are the generators of the fundamental representation of the 
group SU{2). The Hamiltonian depends on the gluon link operators through the "real" and 
"imaginary" parts of the plaquette operators Uu, f/-fc, which are obviously 



^ Tr\U., 



Im (f/i. 



2i 



for 



Uij{x) = Ui{x)Uj{x + ei)Ul{x + ej)U]{x) , i,] = 1, 2, 



(17) 



(18) 



We use Im \Uij\ as an abbreviation for the antihermitean part of the plaquette (which is 
traceless for SU{2)) Utj and Re [t/jj] (which is a multiple of the unit matrix for SU{2)) to 



represent its trace. The Hamiltonian also contains the dimensionless lattice chromo-electric 
field strength operators n"(x), which are canonically conjugate to the links. They obey 
the lattice commutation relations which follow directly from the continuum commutation 
relations in Eq. flTU]) . 

[Utix), U,m = y U,ix) 5,^5,, . (19) 

The constant fj is the product of the near-light-cone parameter rj and the anisotropy param- 
eter ^ = a_/a±, 

v = vi- (20) 

If one chooses rj = \ and varies ^, one simulates an anisotropic equal time theory with a 
ratio ^ = a^/a±_ of lattice constants a_ and a± in longitudinal and transverse directions. In 
the limit ^ ^ one ends up with a system, which is contracted in the longitudinal direction. 



Verlinde and Verlinde 



28 



and Arefeva 



29j have advocated such a lattice to describe high 



energy scattering. A contracted longitudinal system means that even the minimal momenta 
become high in longitudinal direction which is a promising starting point for high energy 
scattering. It is obvious that this limit leads to the same physics as the light cone limit with 
equal lattice constants in longitudinal and transverse directions while ?7 — >■ 0. In both cases 
the near-light-cone Hamiltonian is dominated by the terms proportional to 1/ff involving 
transverse chromo-electric and chromo-magnetic fields. 



In Ref . 



2l| we have determined a variational gluonic ground state wave functional |\I'o) 



which consists of a product of single-plaquette wave functionals with two variational param- 
eters p and 5 



l^o) = *o[f/] |0) = v^e^f^' |0) 



Re U.k(x) 



(5o(A,r7)Tr 



Re U^oix 



(21) 



fc=i 



N^, is a normalization factor. Here, the state |0) represents the trivial ground state which 
is annihilated by the field momenta H^(a;) conjugate to the links, 



H^(f) |0) = and (0| H^(f) = for all x,k,a. 



(22) 



This ground state wave functional is similar to the ground state wave functional used in the 
strong couphng limit of equal-time quantised lattice gauge theory [30|]. However, it takes into 
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account the anisotropy of the gauge dynamics in the purely transversal and the transversal- 
longitudinal planes. As in the equal time case, keeping the wave functional restricted to the 
one-plaquette form does not allow to perform a continuous approach to the continuum limit. 



Further possible improvements are discussed in Ref. [3l|. We have optimized this ansatz 
with respect to the expectation value of the Hamiltonian over a large region in coupling 
space. We are in the position to extrapolate the parameters po, to the light cone fj 0. 
This limit yields the following remaining dependence on A 



Po(A,0)= (o.65-^ + l|^)v/A, 



5o(A,0) = (0.05 + ^-l#^)v/A. (23) 



A A2 

A typical value of A used in the subsequent calculations is A = 10, for which one obtains 

Po(10,0) = 1.83 

5o(10,0) = 0.13. (24) 

Note that the ground state Eq. fl2T|) is an approximation found for the fully interacting 
effective Hamiltonian and does not rely on any truncated Fock space expansion around the 
perturbative vacuum. 

In the light cone limit fj oi the ground state wave functional Eq. (pTI) . under averaging 
with the weight |\E'o[f^]P) the behavior of the gauge fields in the transverse (1,2) plane is 
strongly coupled as shown by the parameter (5o(A,0) becoming small (cf. Eq. (1241) ) . As 
a consequence of this, in the limit (5o(A,0) — > 0, a strong coupling approximation turns 
out to be valid also in the transversal-longitudinal plane even for values of A which are 
far from the region A << 1, say A = 10, as proven by actual Monte Carlo sampling of 
the squared ground state wave functional Eq. fl^T]) [31]]. In the light cone limit fj ^ the 
gauge dynamics in each of the hyperplanes (— , 1) and (— , 2) becomes two dimensional (and 
decoupled). In two dimensions with free boundary conditions, moreover, the strong coupling 
approximation is exact. These are the reasons why Hamiltonian gluon dynamics on the light 
cone is considerably simplified compared with equal time Hamiltonian QCD. We have the 
following standard area law behavior for Wegner- Wilson loops W in the (x~, 1) and (x~, 2) 
planes, namely for 

W{0, 0± ; Z-, di_) = ^_(0, 0± ; z-, 0±) S±{z-, 0± ; z' , d±)S.{z-, di_ ; 0, d±) S^{0, d± ; 0, 0^) 

(25) 
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FIG. 1: In the light cone hmit of the nlc Hamiltonian, the ground state wave functional decou- 
ples the dynamics of the {x~,xi)- and (x^, a;2)-planes from each other. Thus, one obtains many 
decoupled two-dimensional gauge theories. The planar structures (Wegner- Wilson loops) shown 
in black sketch vacuum fluctuations described by the ground state wave functional inside of the 
x-'^)-planes. 

expressible through the plaquettes 

(*o|^Tr[iy(0,Ox;;2-,ci"l)]|^o)= ((^o|^Tr[[/_,]|M/o)^ ' . (26) 

The physical area A of the Wegner- Wilson loop is given hj A = d±\z~\a±a^. Factorization 
is also true for expectation values of the product of two Wegner- Wilson loops which do not 
overlap. Single plaquette expectation values with respect to the ground state wave functional 
are given by 

A. ^ (*o| ^Tr[f/_,]|M/o) = + O(5o') €[-!,!] 



/...(v.o|^^Tr[f/_.]J l^o)^^^^^ + ^ + 0(.,^)e[0,l], (27) 

taking into account that higher powers of the same plaquette do not factorize. Here, 
denote the modified Bessel functions of the first kind. The very same relations hold for 
the purely transversal plaquette expectation values with po substituted by 5o- By using 
Eq. f l26l) and Eq. ( l27j) . one can analytically evaluate all the gluonic matrix elements we need 
for our calculation. One can estimate the physical value of the transversal lattice spacing 
by identifying the rate of the exponential fall-off of a purely transversal Wegner- Wilson loop 
with the dimensionless string tension. For example, one obtains a transversal lattice spacing 
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of a± ^ 0.65 fm at A = 10. This corresponds to a momentum scale of ~ 1 GeV^ which 
is the typical input scale for phenomenological parameterizations of parton distribution 
functions. 

IV. MODELING A COLOR DIPOLE 

The near-light-cone Hamiltonian in Eq. (fT6l) contains only gluon fields, therefore we 
cannot derive hadronic wave functions from this Hamiltonian. We have to make a model for 
the hadron taking into account the gluon structure as exactly as possible while treating the 
quarks only schematically. Our model consists of a dipole state with a fixed longitudinal 
center of mass momentum p_ localized in transversal configuration space at a fixed center of 
mass position 6_l = while the quark and antiquark positions are fixed at ±d±/2, i.e. they 
are separated by the vector d±_ and connected by a Schwinger string along some path C± in 
the transversal plane specifying the dipole in transversal configuration space and longitudinal 
momentum space 



For simplicity, we consider scalar QCD with a scalar matter field. The scalar quark fields 
can be expanded in terms of creation and annihilation operators 



$+(x)= d~ka{k)e-'^'' , = d~kb\k)e+"''' , d~k = j^27r6^^\e - m^) . 



Here, $+(a;) denotes the positive frequency part of the scalar field and ^^{x) represents the 
negative frequency part. The operators a{k),a^{k) refer to the quark annihilation/creation 
operators whereas h{k),h'^{k) refer to the antiquark annihilation/creation operators, m? 
denotes the quark/antiquark mass squared. 

In order to construct such a dipole state, we start with a dipole state which is localized also 
in longitudinal configuration space. Then the dipole state consists of a quark at longitudinal 
position x~ and at transversal position —d^jl and of an antiquark at the same x~ and 
at transversal position d^jl connected by a Schwinger string along the path C_l in the 
transversal plane in order to achieve gauge invariance. The transversal path of n steps 
is parameterized by the intermediate transversal positions (j = l,n — 1) of the links 




(28) 




(29) 



d'^k 
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passed by the path 

C± : -dj2 ^ y^[l] y^[2] ^ . . . ^ - 1] ^ dj2 . (30) 

The entire dipole state locahzed in full configuration space is created by some operator 
acting on the vacuum state 

\Q) = |$o) ® l^o) (31) 

of the entire Hilbert space including gauge fields and (scalar) quarks. The vacuum state |$o) 
of the (heavy) quark sector is assumed to be the Fock vacuum. The vacuum state of the 
gauge fields, however, is given by Eq. fl2Tl) and therefore of non-perturbative nature. The 
operator has the form 

\n)^iix',-dj2)S^^^ix-,-dj2;x-Jj2)[{yT}]c^_(x-,dj2) , (32) 



where the parallel transporter Sg^ represents the path ordered (P-ordered) product of 



n 



transversal link operators along the path C± in the transversal plane. 

In addition, we have allowed in Eq. (l32l) for different longitudinal positions of the transver- 
sal links as motivated below and have inserted longitudinal Schwinger string bits connecting 
the (otherwise adjacent) transversal links in order to retain gauge invariance 



Sflix-, -dj2 ; x~, dj2)[{yT}] = S^x', -dj2- y^, -dj2) S^{y^, -dj2 ; yf, y^[l]) 

' n-l 

'P n ^'(yj^ ^^b1; yj+v y±[j])s±{yj+i, y±[j] ; y±[j + 1]) 



S4y-,dj2; x-,d^/2). 

(33) 



Here, y±[n\ = d±/2. The argument in round brackets of S^g represents the starting and the 
end point of the wiggly string, whereas the argument in square brackets represents the set of 
longitudinal coordinates of the intermediate transversal links. A typical string configuration 
is graphically represented in Fig. [2] assuming a lattice structure. Each of the link operators 
represents a string bit of the entire string. The transversal part of the transporter between 
the quark and the antiquark along the vector d with minimal length (in purely transversal 
direction) represents the ground state of the nlc Hamiltonian in the strong coupling limit. 
In this limit, the nlc Hamiltonian is dominated by the chromo-electric field operators and 
the energy of the dipole state scales with the transversal length \d\ of the gluonic string. 
Because of the Lorentz boost in the longitudinal direction accompanied by the transition 

13 




X +y3 




X 



FIG. 2: Graphical representation of the dipole state on the lattice. The black dots represent the 
quark and the antiquark. The transversal links are allowed to move freely along the longitudinal 
direction. For simplification only transversal links in one direction are shown. The coordinates 
y~ denote the displacement of the link j along the longitudinal direction with respect to the 
longitudinal position of the quark and antiquark. 

from lab frame coordinates to nlc coordinates, the transverse electric field strength operators 
appear in the near-light-cone Hamiltonian with a weight larger by a factor compared to 
the longitudinal field strengths (c.f. Eq. (|T6l) ). Therefore every string with a fixed number 
of transversal links practically has the same energy regardless of the number of links in 
^"-direction in the light cone limit fj 0. This implies that the string should be defined 
having any number of longitudinal links. 

The aim of our calculation is to calculate the gluon structure function for a color dipole 
with string configurations deformed in this manner, however with fixed momentum. We 
have to construct from the localized dipole in configuration space a fast moving momentum 
eigenstate with total longitudinal momentum p„ for which we can determine the momenta of 
gluons. For this purpose we perform an integration over all translations of this state over the 
entire coordinate range multiplied with the appropriate momentum eigenf unctions e"'^-^ . 
Finally, weighting will be performed with the vacuum wave functional squared. Since none 
of the link configurations is preferable from the energetic point of view, we integrate over 
all possible link configurations and assign to each configuration a probability amplitude 
We model by the product of momentum eigenfunctions of each link 
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integrated over all possible link momenta. Due to the projection of the entire dipole state 
onto total momentum p_, the sum of its constituent momenta is restricted (c.f. Eq. f|T^ ). 
i.e. 



^(Ki) = \T\dl'-]e " ^ e U_ - V/M . (34) 




Hence, the final dipole state is given by 

dip^ ; -dj2, C^, dj2)) = II^J] dlidy;^ e" e(^P--J2 (^5) 

X J dx- e-'P-^' ^iix-,-dj2) S^^Hx~,-d±/2; x-Jj2)[{x- + yr}]<^4x-J^^^^ \n) . 

The appropriate normalization of the dipole state is guaranteed by division with a suitable 
normalization factor \/N. This dipole state represents the starting point for our investigation 
of its gluon structure. 

Matrix elements between two dipole states can be computed by contracting the scalar 
operators $± yielding the Feynman propagator Af{x, y\A). We find for the Feynman prop- 
agator of the interacting scalar theory in the eikonal approximation (quark/antiquark have 
large p+ momenta) 

AF(x,y\A) = ^ \e(x- -y~)e-'P-^''~-y'^ + Gfy" - x") 6"^^- 

xS-{x~,x±; y~,x±) 6'^^\x± - y_L) 6{x^ - y^) . (36) 

At high p+-momentum, the quark and antiquark in the color dipole move on straight line 
classical trajectories and pick up non-abelian phase factors along their paths. Instead of the 
usual time ordering, we have an ordering along the longitudinal spatial coordinate. Thus, 
by evaluating matrix elements between two dipole states, additional straight line Schwinger 
strings appear along the longitudinal direction connecting the wiggly Schwinger strings from 
the incoming and outgoing dipole states. To evaluate the expectation value of the point split 
operator G (given in Eq. (ITTl) ) between two dipole states with fixed string configurations 
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the following expression can be reduced to a purely gluonic matrix element 

{n\ (<^l{x'-,x^')S':;!^^{x'-, f/; x'-, x^' + d^')[{y'f}]<i>4x'-,x^' + d^')y 

xG(^^1{x-,x±)S^^{x-,x±; x-,x± + dA_)[{yr}]^_{x-,x± + d±)'^ \n) 
= (^o|Tr S^{x~,x:i + d±; x'~]x± + d±)Sgg^{x'~,x±;x'~,x± + d_i)[{y'f}]^ 
X S-{x'~,x± ; x~;x±)GSgg^{x',x_L ; x~,x_l + d±)[{y~}] |^o) 



to be obtained by averaging over the vacuum wave functional squared. The transversal 

chromo-electric field operators appearing in Eq. ( ITTi) do not commute with the transversal 

link operators appearing in the definition of the dipole operator. Therefore, one has to 

take care of the right arrangement of the operators. The string 5*^^ arising from the 

dipole at x'~ must appear to the left of the operator G{z~ , z± ; z~', zj_) and correspondingly 

the string Sg^{x~) to the right of (^(z",^! ; z~',z±) in the matrix element. The resulting 

Eq. fl37j) allows to express the expectation value of the momentum density operator between 

dipole states by the gluonic vacuum average of the trace over a non-rectangular Wegner- 

c c t 

Wilson loop whose edges are given by the transversal parallel transporters Sg^, S"^^ and 
the longitudinal straight line Schwinger strings connecting the two dipole states. Such 
a Wilson loop is shown in Fig. [3l for simplification with only one transverse dimension 
and the a;~-direction. The full curves represent the strings which connect the quark and 
antiquark in each of the dipole states. The dotted strings arise due to the elimination of 
the quark/antiquark operators in the eikonal approximation of the quark propagator. The 
blue curve corresponds to the point split operator G. Note, that in general the wiggly string 
can extend into both transverse directions and the ^"-direction. After the quark fields are 
eliminated, the x~ integration in Eq. fl35|) can be performed for the incoming and outgoing 
dipole state and yields together with the quark and antiquark momenta from Eq. (!37|) a 
delta distribution setting the sum of transversal link momenta equal to the total string 
momentum 

s(^P--J2i-^ , pi=P--p'L-jt. (38) 

In order to compute the normalization of the dipole state, one simply has to substitute the 
point split operator G in Eq. ( 1371) by the unit operator. 
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FIG. 3: Graphical representation of the generalized Wegner- Wilson loop generated by the two 
color dipole states with a quark q and an antiquark q connected by n-links in transversal direction. 
For simplification only transversal links in one direction are shown. The dash-dotted insertion in 
the so-formed Wegner- Wilson loop represents the gluonic two-point operator with a lattice electric 
field operator at longitudinal coordinate and a lattice magnetic field operator at longitudinal 
coordinate z'~ . 

V. NEAR-LIGHT-CONE GLUON CORRELATION FUNCTION ON THE 
LATTICE 

We come now to the practical evaluation of the lattice counterpart of the gluon distribu- 
tion function. On the lattice, a direct simulation of the average gluon momentum fraction 
becomes subtle. Being the lattice generator of longitudinal translations, the longitudinal mo- 
mentum operator induces translations by any multiple of the lattice unit. We discriminate 
between longitudinal lattice momenta and eigenvalues of the longitudinal lattice momentum 
operator in the following. Similar to the dispersion relation for fermions on the lattice, one 
has for each positive valued eigenvalue in the spectrum of the longitudinal lattice momentum 
operator two possible lattice momenta corresponding to this eigenvalue. Even worse, the 
largest possible lattice momentum corresponds to an eigenvalue of the longitudinal momen- 
tum operator close to zero, far away from its maximal possible value. Hence, by choosing 
the hadron to have the maximal lattice momentum, one finds gluon momentum fractions 
which do not add up to unity (neglecting quark momenta). 

We discuss first how the problem arises and second how to circumvent it. We start 
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discretizing the point split operator of the nlc correlation function Eq. (fTTj) : 



G{z , z± ; z' , z±) 



^lm(^U.kiz'-,z^] 



(39) 



+ 2 Ut{z'-, il) S^,{z'-, il, 2;-, il) Tr 



+ h.c. . 



The field strength in lattice form is expressed here in such a way that Eq. flTT]) follows 
in the naive continuum limit. Here, f/_fc(x) is an average over the "forward" plaquette 
jjtoTw^^^ right of x~ and of the "backward" plaquette U^^^^{x) left of x~ , both — plaquettes 
beginning in x and adjacent to the transversal link Uk{x) : 



with 



and 



U.,{x)^l{U'Tix) + U^.t^ix)) 



(f ) = f/_ (f) f/fe (f + e_ ) ?7l (f + Cfc) (f ) 



f/^,='^-(x) = t/fe(f)t/l(f + Cfc - e.)Ulix - e_)f/_(f - e_) 



(40) 



(41) 



(42) 



Note, that the orientations of the forward and backward plaquettes are the same such that 
the projection of the traceless antihermitean part ImU^k onto a°'/2 becomes proportional 



to F'tk{x) in the continuum limit. 



For z' = z , the point split operator reduces to the dimensionless nlc momentum density 
operator (c.f. Eq. ^) with z = (z~,ij_), the lattice form of which is 



VJz) 



5^i(2n«(f)Tr 



cr 



+ 2Tr 



ylm( U^kiz) 



ylm(f/_fc(£) 



(43) 



which becomes the total longitudinal momentum operator P when summed over the entire 
lattice. The variational ground state wave functional ansatz Eq. (12T!) is an exact eigen state 
of the longitudinal momentum operator with eigen value equal to zero for 5q = 0, i.e. 



(44) 



To find the spectrum of the total longitudinal momentum operator, one has to know how it 
acts on link operators. The commutator of the total momentum operator P_ = X^^'P-(^) 



18 



with a transversal link Uj (y) being part of the gluonic string forming the hadron state gives 
(c.f. Eq. (USD) 



1 [U.{y) U,iy + e_) uUy + e,) - U'_Y\y) UM + f/^' 

-f/l(y - e_) Uj{y- e_) U^y + Cj - e_ 
Xy)U,iy + e^)uUy + e,) 

Uliy- e_) [/,(y - e_) f/_(f + e,- - e_)) + ^(a^) . 



backw / -A 



(45) 



The first line of Eq. fH5l) (containing also "curled-up" plaquette insertions into the gluonic 
string) is exact and will be used in subsequent calculations. The exact commutator of P_ 
with the link Uj symbolized by the arrow t has the following graphical representation 

tUir3 -fD +0-0. (46) 

y + e- y + e- y-e^ y - e _ / 

In the second line of Eq. ( H5l) . we have expanded the result in powers of the lattice spacing up 
to quadratic corrections (9(a^), where stands for (a^, a_ a_L, a^). Then the commutator 
of P_ with the link Uj corresponds to gauge invariant forward and backward translations 
along the longitudinal direction 

pjl4(3 -C 

y ] ^ ^ V y + S- y — 

i.e. leads to a discretized covariant first derivative of the link implemented in a symmetric 
way. The Heisenberg equation of motion for the transversal link on the lattice identifies the 
longitudinal momentum operator as the generator of longitudinal translations. In A_ = 
gauge, i.e. for f/_ = 1 the covariant derivative reduces to an ordinary discretized derivative. 

Eigenstates of the longitudinal momentum operator can be found as sums of transverse 
links along the longitudinal direction modulated by appropriate phase factors up to correc- 
tions quadratic in the lattice constants {a-,a±) 



(47) 



J2V-iz),V,{p.,y^) 



sm{p_)Vj{p.,y±) + 0{a^) 



with 



V,{p.,y^) = J2e-'''-''U,{y-, 



y±) 



(48) 



(49) 
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They arise from projecting transversal links localized in configuration space onto a definite 
longitudinal momentum p_ and are not elements of SU (2) because they are superpositions 
of link operators. The longitudinal lattice momenta must be an integer multiple n of 2 7c/N_ 
with n < < — 1, since the longitudinal light cone momentum for an on shell particle 

is always positive (c.f. Eq. f[T^ ). The momentum p_ of the target is chosen as the largest 
momentum in order to have the maximum resolution in the gluon distribution function 



3d, 



33| 



p_ = ^(iV_/2-l). (50) 



Longitudinal lattice gluon momenta have the resolution 

Ap^ 2 



(51) 



p_ N_-2 

In order to have a high resolution, the extension of the lattice in the longitudinal direction 
has to be very large. 

Eqs. ( H8|50l) imply that the largest lattice momentum p_ yields an eigenvalue of the 
longitudinal momentum operator approximately equal to zero. Even though the gauge field 
is a bosonic degree of freedom, the eigen value sin(p_) of the discretized momentum operator 
looks "fermionic" in the Brillouin zone G [— vr, tt] of the longitudinal momentum, i.e. the 
map p- — > sin(p_) is not injective. 

In order to make it injective and monotonically increasing, we perform a similar but much 
simpler operation as for Kogut Susskind fermions, i.e. we block links on a sublattice with 
half of the lattice spacing along the longitudinal direction. Even sites on the fine lattice can 
be identified with lattice sites on the original lattice. Odd lattice sites on the fine lattice 
lie between two neighbouring original lattice sites. The physical extension Lphys and the 
physical momenta are kept fixed during the transition from the original to the fine lattice 
along the longitudinal direction 

-^phys = N_ a_ , 



2 



1 



a_ = 2a_ , 

We denote quantities on the fine sublattice by superscripts /. Using the fine lattice, we 
define a new momentum eigenstate Vj{p-,y±) on the original lattice by a modulated sum 

20 



over fine lattice links Uj (yj: , y±) 



VAp^, y^) = J2 Uj (yj, y^) ^ . (53) 

Thus, by keeping the physical momenta fixed, the allowed fine lattice momenta are given by 
one half of the original lattice momenta. Thereby we reduce the possible lattice momenta 
on the fine lattice by a factor of two and obtain a one-to-one correspondence between the 
lattice momenta and the eigenvalues on the fine lattice. The right hand side of Eq. (135]) is 
obviously an eigenstate of the longitudinal momentum operator on the fine lattice albeit with 
eigenvalue sin(p'^). However, Eq. ( |53ll is not an eigenstate of the momentum operator on 
the original lattice, because the original momentum operator applied to the block averaged 
state does act solely on fine lattice links at even longitudinal fine lattice sites. Therefore, 
we have to introduce a block averaged longitudinal momentum density V- on the original 
lattice, acting on even and odd fine lattice sites. It is given by the following sum of fine 
momentum density operators VL{zJ , z±) which are defined as in Eq. (1431) with all operators 
on the fine lattice (11^ n^'^ and Uj Uj) 

V.{z-,z^) = 2 (^^V[{2z- -l,z^)+V[{2z-,z^) + ^vL{2z- + l,z^)^ . (54) 

The factor two in front of the definition originates from converting the fine lattice operator 
into an operator on the coarse lattice, i.e. P_ = 2 pI similar to Eq. flS2]) . The effective 
longitudinal momentum density on the original lattice Eq. (|54l) has contributions from even 
and odd sites on the fine lattice. Since we have symmetrized the operator with respect to 
the odd lattice sites on the fine lattice by using one half of the forward and one half of the 
backward contribution, Vj{p-,y±) is an eigenstate of the effective longitudinal momentum 
operator on the original lattice with eigenvalue 

P.,V,{p.,y^)] = 2 sin(p_/2) V,ip.,y^) + 0(aj) . (55) 

Now, the largest possible lattice momentum does also correspond to the largest possible 
eigenvalue of the momentum operator and the eigenvalues are monotonically increasing. 
The above expression Eq. ( 155|) of the longitudinal momentum also appears in the dispersion 



relation for bosons lu = ^/YliC^ sin(pj/2))^ + and defines a one-to-one mapping of the 
lattice momenta Pi in the first Brillouin zone to the energy states. It represents an impor- 
tant stratification which allows to calculate momentum fractions. We define an effective 
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correlation function G{z ,z±;z zj_) on the original lattice by averaging the correlation 
functions G-^^^zJ , z± ; zj', z±) on the fine lattice Eq. (139!) 



z-',z±) = 2 ( ^G^(2^- - 2 z"' - 1, il) + G^(2 2", ; 2z-',z^) 



+ ^Gfi2z- + l,z^;2z-' + l,z,_)^ . 



(56) 



This definition is in agreement with the longitudinal momentum density operator Eq. ([51 
on the coarse lattice for z~ = z~' . Finally, the lattice definition of the gluon distribution 
function is given by 

1 1 



9{P-) 



(57) 



z~,z'~' z± 





,d±/2)^G{z , z± ; 


z' ,z±) 


d{p- ; -d±,C±,d±/2)) 

1 c 


{d{P- 


; -c{L,C_L,f{L/2) 


d{p. ; -d±,C±,d^/2)^ 



In order not automatically to enforce Bjorken scaling, we prefer to express the gluon dis- 
tribution function in terms of the gluon momentum p^_ instead of the momentum fraction 
xb- Here, the gluonic component of the hadronic target state has to be defined by the block 
averaged momentum eigenstates given in Eq. (1531) . 

On the lattice, the following orthogonality relation holds for positive gluon momenta 
p- 



=0 



1. 1 



-11 



(58) 



Taking the real part in the orthogonality relation is sufficient since the point split operator 
is symmetric with respect to an interchange of z~ and z'~. In addition to the continuum 
result 6z~fl/2, there is also a finite size contribution which vanishes like l/N.. One finds 
the average gluon momentum {p^_) 

p- 



Z~ Zj_ 



EE 



(^d{p_ ; -d±,C±,d±/2) 


P-{z- 




d{p_ ; -d±,CL,d±/2)) 

1 c 


(d{p- ; ~d±,C^,d^/2) 


d{p- 


; -c4,Cx,rf"l/2)) 



(59) 



(d{p_ ; -d±,C^,d^/2) 


G{z ,z^ ] z , z^_) 


d{p_ ; -d±,C±,d±/2)) 

1 c 


(d{p_; -d±,C±J±/2) 


d{p- ; -d±,C±,d±/2)'j 



22 



Here, the factor 2ti/N^ is due to the discretized measure of the momentum integration. In 
the infinite volume hmit, the leading contribution is of order 0{1) due to the normalization 
of the dipole state. 

VI. GLUON STRUCTURE FUNCTION OF A ONE-LINK DIPOLE 

We start with the computation of the gluon structure function for the one-link dipole. 
Later we will consider the gluon structure function of the more sophisticated multi-link 
dipole. For a computation we use the dipole state Eq. fl35l) reduced to a single link. Since 
the pure glue Hamiltonian of Eq. (fT6|l does not control quark dynamics we have to choose 
between two alternatives: 

• let the quark and antiquark simply follow the gluon link to which they are attached 
to and fix the quark and antiquark momenta to the correspondent link momentum. 

• impose the quark dynamics of the color dipole externally. Since the total hadron 
longitudinal momentum is given by the sum of the momenta of its constituents, the 
total gluon momentum is then fixed. 

We follow the second alternative and take the mean gluon momentum from experiment. At 
the input scale ~ vr^/a^ = 1.5 GeV^ corresponding to A ~ 10 (c.f. Sec. IIII|) . we use 
the MRST-parameterization (sl and assign a mean momentum fraction = 0.38 p_ to the 
string. The string momentum results from the difference of hadron momentum and quark 
and the antiquark momenta which is taken from experiment: 

jP_ = "P- — -pL — jp_ . (60) 

We ascribe this momentum p_ to the complete string of gluon links. Its transverse size 
now equals one of the lattice unit vector Cj, where j = 1, 2 denote the transversal directions 

c4 = e,- , J = 1,2. (61) 

According to Eq. flTTI) . the norm of the one-link dipole is related to the matrix element of a 
Wegner- Wilson loop defined by the eikonal trajectories SL-, SL^ of the quark and antiquark 
(dotted lines) together with the strings (full lines) S^q^\x'^\^ ^q^{^~j \ inside of the color dipole 
(c.f. Fig. [3]). For a string consisting of a single link, we can not define the dipole state 
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symmetrically with respect to the origin in transversal space. We choose without loss of 
generality the quark to be located at the origin in transversal space and extend the dipole 
along one of the positive transversal axes with the antiquark located at ej. Hence, the norm 
of the dipole state is given by 

X Sg{x'f, 0^ ; x'f, dj_) [x'f] S[{x'f, Ox ; x] , 0^) S!^^{x] , 0^ ; xj , dj_) [xj]\ |^o) -(62) 

Here, we have suppressed the path index C± in the state vectors. We absorb factors 2 7r, 
volume factors due to the squared delta distributions and the factor p'Lp'L/m'^ appearing in 
Eq. (l37Il into the common normalization factor A^. In the strong coupling approximation, 
this reduces to 



d{p_ ■ Ox, d^) d{p. ■ Ox, > = ^ Hp-), 



(63) 



where Fi{p^) is given by 



Flip: 



J2 e-'^'-/2(-7-4-)/l 



-Tr 



\X J, X j: I 



(64) 



In the gluon correlation function one has to take care of the arrangement of the opera- 
tors. The transversal chromo-electric field operators in the point split operator Eq. ( l39i) do 
not commute with the transversal link operators appearing in the definition of the dipole 
operator. As in the previous Sec. |Vl the string Sqq[x'j^] arising from the dipole at 
must appear to the left of the point split operator G{z~ , z± ; z~', z±) and correspondingly 
the string Sgg[xj] to the right of G{z~,z±_ ; z~', z±) (see Fig. [3]). Then the forward matrix 
element of G{z~ , z± ; z~', z±) is given by 

d{p- ; Ox, d^) G{z', z±_ ; z~', za_) d{p^ ; Ox, d^)) 



^1 -ip?/2(x7-.-) 

,Xj 

(^ol SL{xJ,d±; x'f;d±) S^j{x'f,6±; x'f , dA_)[x'f] SL{x'~ ; x~,6±) 

G{z~,z±] z~',z±) S^q{xj,0±] xj,d±)[xj] ^ l^o) ■ (65) 

The square brackets denote matrix elements with color indices ab and ba respectively such 
that the expectation value is given by the trace over the product of the color dipole states 



24 



Sqq and Sqq with the effective point split operator G in between. The momentum correlation 
function Eq. fl65l) evaluates the cross product of electric and magnetic field strengths sepa- 
rated along the light cone, i.e. it determines the correlation of an electric field in the dipole 
with the corresponding magnetic field. In order to compute it, we arrange the operator G 
with |\E'o) (c.f. Eq. (1211) ) in a way such that G stands directly in front of the trivial ground 
state |0) 



\0) + S;,^oG |0) . 



(66) 



The trivial ground state is annihilated by this operator. The commutator of G with the 
ground state wave functional leaves the dipole operator intact and yields a vacuum transition 
which is subtracted when the connected matrix element is extracted. Therefore, the only 
remaining contribution comes from the commutator of G with the transversal link of the 
incoming dipole. It is given by 

{zJ ,zj_; zj',z_i_) , Uj {xJ , f _l) = 

^Si{xJ,x±; zJ ,x±)lm(^Uij{zJ ,x_l)^ Si{zJ,x±; xJ,x_L)Uj {xJ ,x±) 6^-^^-, 6g^^g^ 



+ [Zj ^ Zj 



(67) 



Due to the interchange symmetry z ^ z' of the commutator, only the cos-part of the 
Fourier transformation survives: 



N-/2-1 



E 



~ip^_{z —z' ) 



z-,z-'=~N-/2 
N-/2-1 



-N-/2 



G{z ,z^\ z ',zj_) , U^{xj:X_l) 

2 cos (p^ z~) SL{x~ ,x±; xJ + 2 z~ ,x±)lra(^U[j{xJ + 2 z~ ,x±)^ 

X SL{xJ + 2z~,x±; xJ,x±)Ul{xJ,x±)5s^^g^. (68) 



The gluon distribution function Eq. (1581) for a one-link dipole with total string momentum 



becomes 



(^o|Tr [5^(0,rfl; x7,dl)f//^(a;7,0^)5^(x7,0x; 2^-,0^) 

xIm(f/4.(2;2-,0x))s^(2;2-,0x; 0,0x)f//(0,0x)] |^o) • (69) 



25 



In gi{p^_^;p^) we indicate the total string momentum by the argument and the number 
of transversal links by the index n = 1. 

In the following, we discuss the evaluation of the matrix element in Eq. fl69l) . The imagi- 
nary part of a Uij plaquette (field strength) located at longitudinal position 2 has to be 
correlated with a closed loop of links in the longitudinal transversal plane located between 
longitudinal positions and xj as visualized in Fig. HI The field strength is connected with 



2z- 



f — >—• — > • > — • • — • I — • — ► • ► 
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FIG. 4: Visualization of the matrix element in Eq. (j69p . There are two different possible situa- 
tions resulting in different expectation values of the operator, namely either the field strength lies 
outside of the Wegner- Wilson loop (a) or it lies inside (b). The difference in each of the two cases 
corresponds to the antihermitean part of the plaquette. 

the edge of the Wegner- Wilson loop at longitudinal position 0. We distinguish two cases 
depicted in Fig. H] which can be fully evaluated in the strong coupling approximation: 

(a) The plaquette U^j lies outside of the loop. The matrix element factorizes and vanishes 
due to TT[lm{Uij)] = 0. 

(b) The plaquette lies inside the loop, the matrix element can be computed by tiling which 
yields the second line of Eq. ( ITOl) . 



One finally obtains for the gluon distribution function of a one-link dipole 

z 



TT / COS (P_ Z ^ , „, > Sm Xf\ 



1 



111. 

2 



ill. 

2 



\x-\-l 



X 



]^{Qq{2 z-)Q{x-j - 2 Z-) + Q{2 z-)Qq{x-j - 2 z 
-]^{q{-2 z-)Qo{2 z- - x-f) + Qq{-2 z-)Q{2 z 



■ (70) 
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FIG. 5: With crosses we show the expectation value of the longitudinal momentum operator P_ of 
a one-link dipole state as a function of p- = on a = 100 lattice for A = 10 and 100. We also 
show the expected eigenvalue 2 sm{p^ = p^/2) of the longitudinal momentum operator without 
O(a^) corrections of a single- link color dipole state projected onto longitudinal momentum p_. 

The two cases a) and b) are encoded in the G-distributions Oo(a;) and B(x) differing in the 
value they take at x = 0, i.e. 



Let us now check the gluon momentum sum rule by evaluating the expectation value 
of the longitudinal momentum operator (pi) = Xl^- ^ as a function of the 

external string momentum p'^. In Fig. Owe show the mean gluon momentum as a function 
of the lattice momentum p^. The expectation value is obtained on a = 100 lattice 
for two different values of A, i.e. A = 10 and 100. In the region of small momenta p^ 
one can recognise that the lattice discretisation is very accurate. Due to the introduction 
of the finer sublattice the mapping of the lattice momentum p^ to the mean momentum 




(71) 



In the first line of Eq. ( ITOl) . only the imaginary part sm{p^/2xj-) survives from the 
exponential e"'^-*^^^/ due to the antisymmetry of the sum in xj . 
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is unique. For comparison we also show the exact eigenvalue 2sin(p;^/2) without order 
O{o?) corrections. The larger the lattice coupling constant, the more accurate becomes the 
mapping from lattice momenta to observed momenta. 

The full gluon distribution function p^__ gi{p^__^; p^) of a color dipole with one link in the 
transversal direction is shown in Fig. O It is computed on lattices with N_ = 20, 30, 50, 
and 100 sites in the longitudinal direction. The average gluon momentum (p^) of the dipole 
state has been adjusted to the average gluon momentum (p^ ) = 0.38p_ of the MRST gluon 
distribution function at = 1.5 GeV^. So far, the simulated gluon distribution function on 
the lattice depends on the total hadron momentum. This is in contrast to Feynman scaling, 
where the gluon distribution function is only a function of the gluon momentum fraction 
P-/p-- For fixed A, the hadronic lattice momentum p_ = 2tt{N_/2 — 1)/N_ is exclusively 
determined by the longitudinal lattice extension A^_. Independence of the gluon distribution 
function on the hadronic momentum would be equivalent to independence on the lattice 
extension. This figure demonstrates the effect of increasing the number of longitudinal 
lattice sites, i.e. approaching the infinite volume limit. Scaling for A = 10 seems to be 
obeyed for longitudinal lattice extensions larger than = 50. Realistic lattice simulations 
with an improved ground state wave functional need quite large longitudinal lattice sizes. 
The smearing of the distribution function is due to the gluon dynamics incorporated in the 
Wegner- Wilson loop expectation value. The area law behavior of the Wegner- Wilson loop 
yields a non-trivial gluon wave function which broadens the distribution. 

If one varies the lattice gauge coupling A, the single plaquette expectation values vary 
between and 1 depending on the coupling constant A. This has consequences for the 
width of the gluon distribution function p^_ gi{p^-', p'l) as shown in Fig. [71 The larger A, i.e. 
the smaller the QCD gauge coupling g"^, the stronger the peak in the one-link distribution 
function becomes. 

In the extreme weak coupling limit A ^ oo the single plaquette expectation value ap- 
proaches 1. For finite coupling constants, the single plaquette expectation value is less than 
one. Hence, it suppresses large Wegner- Wilson loop extensions in Eq. fl7U]l : 




(72) 
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FIG. 6: Gluon distribution function ; P- = 0.38 p_) of a color dipole with a single transver- 

sal link (n = 1) in the transversal direction for different lattice sizes N^. The distribution functions 
are computed with a lattice coupling A = 4/^^ = 10 in the effective SU{2) lattice Hamiltonian 
of Eq. (jl6p . The average gluon momentum {p^__) = of the dipole state has been adjusted 
to the average gluon momentum {xb)p- = 0.38 p_ of the MRST gluon distribution function at 
Q2 = l.5GeV2. 

One can define a correlation length which represents the longitudinal distance at which the 
single plaquette expectation value reduces to one half of its original value 

Ae = log(^)/log(^(^Tr 

By using Eq. (1251) . Eq. (ITTI) and Eq. (1751) . we can evaluate the correlation length directly as 
a function of A. At A = 10, the correlation length is given by A^I^^^^^q = 3.12 and at A = 50 
it is given by Ai^I^^^q = 8.05. Since the nlc gluon correlation function only has support 
when it lies inside the Wegner- Wilson loop, the correlation length A^ gives an estimate for 
the width Ap^ of the gluon momentum distribution, 

A/_ = • (74) 



U 



-k 



(73) 
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FIG. 7: Gluon distribution function gi{P-]P^ = 0.38 p_) of a color dipole with a single transver- 
sal link (n = 1) in the transversal direction at different lattice gauge couplings A. The distribution 
functions have been computed on a lattice with = 1000 sites in the longitudinal direction. 
The average gluon momentum {p^_) of the dipole state has been adjusted to the average gluon 
momentum {xb)p~- = 0.38 p_ of the MRST gluon distribution function at = 1.5 GeV^ . 

This implies that the width of the gluon momentum distribution at A = 10 is given by 
^P-Ia=io ~ ^■^'^ and at A = 50 by /\p^_\^^^Q = 0.13 as seen in Fig. [3 In the extreme weak 
coupling limit, when the link reduces to a single gluon, the gluon distribution function is 
sharp, i.e. gi{p^_;p^) = 6{p^_—p^). Since the link momentum is fixed by the projection onto a 
definite momentum, also the gluon momentum is fixed in this limit and no variance is allowed. 
On the other hand, for smaller values of A = ^, i.e. for strong coupling the correlation length 
becomes smaller which implies that one has a broad momentum distribution peaked around 
p^. 

VII. THE GLUON DISTRIBUTION FUNCTION OF A HADRON 

Up to now, we have considered the gluon distribution function of a color dipole consisting 
of a single transversal link. The one-link dipole gluon distribution is the basic building 
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block from which the multiple link dipole gluon distribution function of a hadron can be 
constructed. We expand a hadronic state 



h{p-,0±)j in dipole components, i.e. 
-^h{Cj±) = (d{p_; -dj2,C^,d^/2) h{p_,6^)). (75) 



The hadron has the same momentum p_ and the same transversal cm coordinate x± = 
as the dipoles. The wave function '^h{,C±, d±) represents the probability amphtude to find a 
dipole with a quark and an antiquark separated by the transversal distance d± and connected 
by the path C± in the hadron. Hence, the actual hadronic gluon distribution function arises 
from a superposition of multiple link configurations. The wiggly strings Sgg c.f. Eq. 
connecting the quark/antiquark are not restricted to lie along one of the coordinate axes 
(c.f. Fig. [2]), but have a fixed number of transversal links as explained in Sec. IIV[ In order 
to project this state on angular momentum = 0, we rotate the hadron in the transversal 
plane by summing over randomly chosen curves C± which can be constructed in the following 
way. A random walker starts at an initial time t = and trails a Schwinger string along 
its path through the transversal lattice. In each time step, the walker may hop with equal 
weight in one of the four transverse directions. The random walk ends if the number of 
hops corresponds to the number of allowed transversal links fixed by the energy constraint. 
The starting point of the random walker has to be chosen a posteriori in such a way that 
the center of mass of the generated dipole configuration is at the origin. The ensemble of 
possible random paths automatically obeys the desired rotational symmetry. 

Since the energy of the strings with a given number of transversal links is the same for 
all the string configurations, we assume that the probabilities among the total number 7^ of 
curves with n-links are equally distributed: 

\^h{CJ)\' = ^. (76) 

From the random walk follows that for n-links the hadron has an average radius squared 
proportional to n: Hence, the area of the hadron scales with the number of links 

<^i> = ^ • (77) 
Due to the strong coupling approximation non vanishing gluonic matrix elements need 
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incoming and outgoing states to have the same curve connecting the quark and the antiquark: 

J%_,0^)|o|%_,Ox) 
J2 ; -d^/2,C±,dj2)\0 |%_ ; -dj2,C^,dj2)'j . (78) 

C,d± 

Because of the equal weight of all the dipole configurations with fixed transversal length 
n, the gluon distribution can be calculated from the distribution function gn of a string 
elongated along only one of the transversal axes (c.f. Fig. [2]): 

gh{p'-;p'-)=9n{p1;p'-). (79) 

Due to the sum rule (c.f. Eq. fl59|) ) the expectation value of the gluon momentum inside 
the n-link dipole is fixed as {p^_) = p^. The computation of Qn is done in analogy to 
the computation of the single-link gluon distribution. Because of the summation of the 
nlc correlation function over the entire transversal lattice, the chromo-electric field 11^ in 
the operator defining the nlc correlation function, can act on each of the transversal links 
appearing in the connector Sqg. In strong coupling the total loop factorizes, therefore the n- 
link distribution function is given by the product of a splitting function P„^„_i multiplying 
the gluon distribution function with n — 1 links . In this recursion relation (c.f. app. E]) all 
possible intermediate momenta of the substring are summed over: 

9n{pt;p'_) = ^ E 9n-l{p'-;p'J) Pn^n-liptp'l') • (80) 

p£'=o 

The splitting function Pn^n-i{P-,P-') denotes the probability that a string with n 
transversal links and total momentum splits into a string with n — 1 transversal links and 
total momentum p'^'. The form of the splitting function is of a kinematical and a dynamical 
origin which are both encoded in the functions F^ip^) for m = n — 1 and m = 1 due to the 
recursive representation: 



n 



Fn-l{p'-')F^{p^^-p'_') 



Pn^n-liplp'J) = ( ) X J^-'-^V^-^^-V^- . (81) 

The functions Fn-i{p^) and Fi{p^) contain momentum conservation in the density matrix 

Pn, 

n _■ ([3 12 ^ —l-^ ' l2 ' 

Pn[ptyj. y'n = E ^(P- - E ^(P- - E ^- ) ^" ^ '~ ^ , (82) 
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and the gluon dynamics in the residual part of Fn{p^) is related to the n-fold product of 
expectation values of Wilson loops with one link in transverse direction and \y^ — y'"] links 
in longitudinal direction (cf. Eq. 



Fn{p'.)= E Pn[P-;%>%in(^Tr[t/_.])'^^""^^"''. (83) 
{y-},{y-'} ■'=1 

The denominator of Eq. flSTl) guarantees the correct normalization of the splitting function 
Pn^n-i{P-,P-') which has to satisfy the following relation 

^j2p-'p----^^p-y-')=p-^ (84) 

p1'=0 

in order that the distribution function gn{P-,P-) obeys the momentum sum rule 

^f^f-9n{f-;p'.)=p'.. (85) 



The initial condition for the recursion relation Eq. flHOl) is given the one-link dipole func- 
tion gi{p^^',p'l) derived in Eq. (!70|) with the total gluon momentum fraction taken from 
experiment. We use as lattice gauge coupling X = A/g^ = 10 unless otherwise noted which 
corresponds to typically "strong coupling " transverse lattice size far from the continuum 
a± = 0.5 — 0.6 fm, i.e. to an input scale of^^ ~ Ti^/a]_ = 1.5 GeV^. One can try to devolve 
the phenomenological NLO MRST 2002 [5] and the CTEQ GAB parameterizations of 
the gluon distribution function and one finds p^ = 0.38 Since the computation is purely 
arithmetic in strong coupling , we can use a large longitudinal lattice with = 1000 lat- 
tice sites. The so defined lattice gluon distribution function depends on the gluon lattice 
momenta p^_ = 2 irn/N^ where is the number of lattice sites in longitudinal direction on 
the coarse lattice and the integer < n < N_/2 — l. With = 1000 we find a smooth limit 
for the structure function, which we can associate naively as a scaling structure function 
(cf. Fig. [6]). The simple vacuum wave functional we use does not allow us to discuss the 
continuum light cone limit with the longitudinal lattice size N_ a_ constant when oo 



M 



35| that the infinite volume 



and a_ 0. It has been shown in the Schwinger model 
limit has to be performed before the light cone limit. 

If one increases the number of transversal links, the gluons have access to a larger region 
in phase space due to the splitting function Pn^n-i in Eq- (IHOl) . An increase in the number of 
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transversal link operators implies that the total gluon momentum will be partitioned among 
more gluons. Hence, it becomes more likely to find a gluon with a small fraction of the total 
momentum. This can be observed in Fig. [HI The mean momentum fraction of the gluons 
i.e. the integral under the curve remains constant, however, the gluons with large momenta 
are shifted from large to smaller values of . 
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FIG. 8: Gluon distribution function xb 9nixB]P-) of a color dipole with different number of links 
(n = 1,2,4) in the transversal direction. The distribution functions have been computed on a 
lattice with = 1000 sites in the longitudinal direction at a lattice gauge coupling A = 10. The 
average gluon momentum of the dipole state has been adjusted to the average gluon momentum 
{xb) = 0.38 of the MRST gluon distribution function at = 1.5 GeV^. 



We can discuss what happens when one increases the resolution. Keeping the transversal 
extension of the dipole fixed we increase the number of transversal links by one unit, then the 
recursion relation Eq. flHOl) gives a strong coupling equivalent to the weak coupling DGLAP 
equation, which describes the change of the parton distribution function under a variation 
of resolution. Indeed, if one has scaling in the limit p'^ 00 the recursion relation Eq. flHOj) 
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can be written as 



gniXs) 


— [ gn-liXs/ Zb) Pn-*n~l{zB) , 
JxB 




Xb 


= P-/P-^ 




Zb 


= P^'/p^ ■ 


(86) 



Subtracting gn-i from gn one arrives at an equation which has almost the form of the weak 
couphng DGLAP equation. The main difference occurs in the redefined sphtting function 
Pn^n-i- In the usual DGLAP equation, the splitting function denotes the probability for 
a gluon to split into two gluons, one of them carrying the momentum fraction Zb- Our 
equation resembles more the LUND model |36i where the dynamics of the entire fragmenting 
string is described. One can see that in the weak coupling limit the plaquette expectation 
values become unity plus 0{g^) corrections and make the redefined splitting function P — 1 
proportional to a^. Once the continuum limit is under control with a suitable wave functional 
of the ground state, one may consider the transition of the so redefined splitting function 
into the DGLAP Kernel. 

In Fig. [9|, we compare the theoretical gluon structure function for a n = 4 link dipole with 
the MRST and the CTEQ gluon distribution function at = 1.5 GeV^ as functions of the 
gluon fractional momentum xb = P-/p-- As before, the first moment of the lattice gluon 
distribution function has been fixed in this figure to the value (xb) = 0.38 at = 1.5 GeV^. 
The average gluon fractional momentum obtained from the CTEQ parameterization differs 
only by ten per cent from the MRST value. We choose four links to be consistent with the 

2 

size of the proton and the relation {R'j_) = and a transversal lattice size of a_L ~ 0.65 fm. 
The functional behavior of the gluon distribution function as a function of multiplied with 
is the same as the functional behavior of the gluon distribution function as a function of 
■p^_ multiplied with pi , 

XB9n{xB\P^-) = p'Lgn{P-;P^)\pa_^^^p_ ■ (87) 

The lattice gluon distribution function agrees within the systematic uncertainty with the 
phenomenological MRST -gluon distribution function . But the figure shows that there is 
a large systematic uncertainty in the gluon distribution function evolved to = L5GeV^ 
depending on the different parameterizations. The MRST collaboration even gives negative 
values of the gluon distribution function at small values of x^ for such a small Q^. 
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FIG. 9: Gluon distribution function xb gnixB',P-) of a color dipole whose number of links in 
the transversal direction is given by n = 4 in comparison with the MRST and the CTEQ gluon 
distribution function at = 1.5 GeV^. The lattice distribution function has been computed on a 
lattice with = 1000 sites in the longitudinal direction at a lattice gauge coupling A = 10. The 
average gluon momentum of the dipole state has been adjusted to the average gluon momentum 
{xb) = 0.38 of the MRST gluon distribution function at = 1.5 GeV^. 

An important property of the gluon distribution function at small values of is its 
dependence on hadronic size. One knows hadronic cross sections at intermediate energies 
and can deduce that the gluon structure function at small or the soft Pomeron coupling 



depends on the area of the hadron 



371]. With decreasing xb the gluons become uniformly 



distributed inside the hadron such that the gluon distribution function should indeed be 
proportional to the transversal area vr(i?^) of the hadron. 

In Fig. [T0|, we show the gluon distribution function at the lowest value of xb compatible 
with the lattice momentum cut-off, i.e. xb = = 0.002 as a function of n for two different 
values of A: 



X 



B 



—r 

For A = 50, the gluon distribution function at x™"^ depends linearly on the hadronic size 
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n — 2 i.e. one obtains the expected dependence of the "hadron cross-section" . In 

order to guide the eye, we also plot the best fit with 

5n(xr ;P-)L=50 = c^ ' ^ = 0-42 (89) 

into the plot. 




n 



FIG. 10: Gluon distribution function a;™° 5f„(x^''^;p:^) as a function of the number of transversal 
links n for two different values of A. The average gluon momentum of the dipole state has been 
adjusted to the average gluon momentum [xb) = 0.38 of the MRST gluon distribution function at 
Q2 = 1.5 GeV^. 

For A = 10, i.e. for stronger coupling , the dependence of the gluon distribution function 
at a;^"^ on n is less then a linear. In the strong coupling regime, rotational invariance on the 
lattice is broken. Therefore, the cross section of the hadron is no longer given by a circular 
disk. 

VIII. SUMMARY AND OUTLOOK 

In high energy scattering partons move along almost light like trajectories. Hence, light 
cone coordinates define the appropriate framework in order to describe high energy scatter- 
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ing experiments. If one wants to apply the computational methods of lattice gauge theory 
usually defined in a Euclidean signature to the computation of observables of high energy 
scattering experiments one has to face the problem that the light cone shrinks to a single 
point. Hence, correlation functions along the light cone which are important for the de- 
termination of structure functions can not be computed directly. One needs the operator 
product expansion in order to compute these correlations on the lattice. By doing so, one 
is restricted to the moments of the structure function. 

We have proposed to use the nlc lattice formulation in order to compute the correlation 
functions on the light cone directly. We have generalized the definition of the light cone 
correlation function to nlc coordinates, such that in the light cone limit the original definition 
is recovered and that the nlc correlation function obeys momentum conservation. In our 
approach, we are not restricted to the moments of the gluon distribution function. 

We employ the nlc ground state wave functional in the light cone limit which was varia- 
tionally optimized close to the light cone limit (cf. Ref. 211]). Since our theory is formulated 
in a Hamiltonian framework, we stay in Minkowski space-time throughout the computation. 
This implies that one does not need to perform an analytical continuation from a Euclidean 
to a Minkowskian signature at the end of the computation. 

The nlc ground state wave functional ansatz in the light cone limit shows a significant 
simplification for the computation of gluonic matrix elements in comparison to an equiv- 
alent equal time quantised computation. The ground state wave functional decouples the 
purely transversal dynamics. Hence, one effectively deals with two two-dimensional gauge 
theories, each living in a set of longitudinal-transversal planes which are distinguished by 
the other transversal coordinate. If one neglects boundary terms, the analytical tools for the 
computation of matrix elements valid in the strong coupling approximation become exact 
over the entire coupling regime. 

We insert a color dipole state into the vacuum described by the ground state wave func- 
tional. The construction of the hadronic dipole is guided by the principles of the strong 
coupling approximation, i.e. the Schwinger string connecting the quark and the antiquark 
is chosen to follow the minimal transversal path in between the quark and the antiquark. 
The lattice naturally defines the hadronic state in configuration space. This is in contrast to 



other light cone lattice approaches like the transverse lattice approach 



18| where so-called 



"fat links" are quantised canonically and as such have an explicit formulation in the mo- 
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mentum representation which is more natural for the computation of structure functions. 
In our approach, we project the configuration space states exphcitly on states with defi- 
nite momenta in such a way, that each of the hnks forming the Schwinger string has its 
own momentum. Only the total momentum of the Schwinger string is constrained by the 
total hadron momentum (minus the quark and antiquark momenta). Since our quarks are 
not dynamical, we cannot obtain their momenta from within our calculation. Therefore we 
take them from experiment. We use the average total string momentum obtained from the 
MRST (2002) NLO parameterization of the gluon distribution function at = 1.5 GeV^ 
as an input. 

The so obtained gluon distribution function obeys a recursive equation which relates the 
gluon distribution function with n transversal links to the gluon distribution function with 
n — 1 transversal links. If one interprets the increase of link constituents at a fixed size 
of the dipole as an increase in the resolution of the probe, this recursion relation is the 
non-perturbative counterpart of the DGLAP equation. Indeed, with increasing number of 
transversal links, the gluon distribution function grows at small xb due to the fact that 
the available total gluonic momentum has to be distributed among more and more con- 
stituents. A string-splitting function represents the probability to find a string containing 
n — 1 transversal links inside of a string with n transversal links. 

Our results calculated for the QCD-coupling X = 10 = 4/ g"^ roughly correlate with a 
transverse lattice spacing of a± = 0.5 — 0.7fm. They can be compared with phenomenological 
parton distributions, if we choose the number of links appropriately for the proton (n = 4). 
The calculated low x gluon structure function shows a behavior similar to the MRST- 
parametrization, once we fix the mean (xb) in accord with these. Unfortunately, due to the 
lack of quark dynamics the mean gluon momentum (xb) itself is out of reach. 

The model presented here also shows that xb g{xB) for the gluon at small x becomes pro- 
portional to the hadronic size R\. This coincides with the empirical soft Pomeron behavior 
of hadronic cross sections. Both, the evolution of the structure function with increasing 
resolution and/or with decreasing xb need a more sophisticated ground state wave func- 
tional (respecting scaling with the lattice spacing) and numerical simulations (corresponding 
to the inevitably non- locale action). In previous work an improved wave functional has been 



proposed 



311 ] which can be used for structure function calculations, once it has passed the 



scaling tests in the light cone limit r] 0. 
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APPENDIX A: RELATION OF TO Qn-i 



The dipole matrix element is related to Fn{p^) which is given by 



With this definition, we first prove the following relation for Fn{p^): 



(Al) 



(A2) 



pt'=0 



One can use the definition of the density matrix and can insert a unity in form of two 
additional momentum summations with appropriate Kronecker deltas in order to obtain 



I \ H r " 1 



IS/, -y, 



n n _ 

j=2 j=2 



E 

i E {l'-/2y~-fj/2y'^-) 



■ (A3) 



The first square bracket in the above equation is evaluated easily by performing a variable 
transformation and is given by 



yi .s/r 



Is/i I „-ip_/2y- 



(A4) 



If one inserts this result into Eq. (1A3I) . one can perform the p'_ summation. After the 
evaluation of S{p- — p'_), one can identify the second square bracket with — p^') 

and Eq. ( ]A2I) is proven. 

By using the same method to split of the contribution of a single link from the entire 
matrix element, one can show that the following relation holds for the gluon distribution 
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function of a dipole containing n transversal links: 

The factor of n in front of the sum is due to the fact that the correlation function successively 
acts on each of the transversal links assembling the n-link dipole state when summed over 
the entire transversal lattice. The factor 1/Fn{p^) ensures the correct normalization of the 
n-link dipole state. The sum goes over Fi{p'^') gi{p^_; p^') times Fn-i{p^ —p^'). 

By using Eq. flA5l) . one can rewrite the recursion relation in terms of Qn-i in order to 
obtain 

9n{p'-;p^l) = (^) E 9n-i{p'-;p'-')F^^,{p'_')F,{p'_-p'_'). (A6) 
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